Complete WMAP Constraints on Bandlimited Inflationary Features 
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Using a principal component (PC) basis that accommodates order unity features in the slow roll 
parameters as fine as 1/10 of a decade across more than 2 decades of the inflationary expansion, 
we test slow roll and single field inflation with the WMAP7 data. Detection of any non-zero 
component would represent a violation of ordinary slow roll and indicate a feature in the inflaton 
potential or sound speed. Although one component shows a deviation at the 98% CL, it cannot 
be considered statistically significant given the 20 components tested. The maximum likelihood 
PC parameters only improves 2AlnL by 17 for the 20 parameters associated with known glitches 
in the WMAP power spectrum at multipoles £ < 60. We make model-independent predictions 
for the matching glitches in the polarization spectrum that would test their inflationary origin. 
This complete analysis for bandlimited features in the source function of generalized slow roll can 
be used to constrain parameters of specific models of the inflaton potential without requiring a 
separate likelihood analysis for each choice. We illustrate its use by placing bounds on the height 
and width of a step-like feature in the potential proposed to explain the glitch at 20 < I < 40. Even 
allowing for the presence of features in the temperature spectrum, single field inflation makes sharp 
falsifiable predictions for the acoustic peaks in the polarization whose violation would require extra 
degrees of freedom. 



I. INTRODUCTION 

Observed glitches in the WMAP power spectrum of 
cosmic microwave background (CMB) temperature fluc- 
tuations [l[ have motivated many studies of features in 
the initial conditions. Most of these studies have focused 
on the reconstruction of the curvature power spectrum 
through parametric, minimally parametric or regularized 
inverse techniques (e.g. 0412 )• 

Reconstruction of the curvature spectrum suffers from 
two potential issues. Given that fine scale features are 
observable at high wavenumbcr, parametric models are 
not complete unless a very large number of parameters 
are employed. Secondly, not all curvature power spectra 
can arise from physical mechanisms in the early universe 
making parametric models potentially overcomplete and 
subject to fitting the noise instead of fitting the physics. 
For example, a delta function in the initial curvature 
spectrum would be highly observable but not expected 
to arise in any physical model. 

For the purposes of testing inflationary models of the 
initial conditions one can instead try to constrain the 
shape of the inflaton potential under the assumption that 
inflation arises from a single scalar field with a canonical 
kinetic term. Specific potentials have been used to test 
the origin of the low quadrupole moment, the glitches 
at multipole moments £ = 20 — 40 and glitches near 
the WMAP beam scale Q, [T|Hl5j]. On the other hand, 
model-independent reconstruction approaches have im- 
plicitly or explicitly assumed a slowly varying inflaton 
potential [16l - F20j . 

Sharp features in the inflaton potential cause features 
in the temperature power spectrum [2l|, [22[ . As long as 
those features are of small amplitude, inflation contin- 
ues uninterrupted but certain slow roll parameters arc 



neither constant nor necessarily small. The generalized 
slow roll (GSR) approach [23l - [26j can be used to analyze 
such cases. In particular to good approximation there is 
a single source function that encodes observable features 
in the inflation potential [25j | for canonical kinetic terms 
or the sound speed for non-canonical terms [27j |. This 
function is also closely related to the source of corre- 
sponding bispectrum features (28|. In previous work, we 
studied the strong constraints on this function imposed 
by the precise and featureless measurements around the 
first acoustic peak through a low order principal compo- 
nents decomposition [29| . 

In this paper, we extend our previous analysis to a 
basis of 20 principal components for the source function 
of inflationary features. This basis is complete for mod- 
els where the features vary no more rapidly than 10 per 
decade of the expansion or about 4 per efold during in- 
flation. In SjlT]we review the GSR and principal compo- 
nents technique. In Appendix \K\ we describe numerical 
techniques used to reduce the computation time of the 
analysis. We test the validity of the GSR approximation 
in Appendix [Bj We present the results of the WMAP 
likelihood analysis in mill In iflVI we develop tests of 
single field inflation and consider applications to specific 
classes of potentials. We discuss these results in SjV] 



II. METHODOLOGY 

We use the generalized slow roll (GSR) approximation 
and principal components to study features generated by 
single-field inflation. We refer the reader to Ref. [29[ for 
details but provide a brief description here. Features arise 
from a single source function that describes the deviation 
from slow roll in the background ( mi A[) which we can be 
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decomposed into a basis of principal components that is 
complete for bandlimited models sampled at a maximal 
rate per efold ( §11 Bj) . 



A. Generalized Slow Roll 

Under the GSR approximation, features in the curva- 
ture power spectrum are generated by a single source 
function of the background evolution of the inflaton </> 

G'Qnr?) ee -2(ln/)' + H(ln/)", (l) 



where 



/ = 



2-K(f>ar] 
H 



(2) 



Primes here and below denote derivatives with respect to 
In 77 where 77 = J t *° nd dt' /a(t') is the conformal time to the 
end of inflation and we take units of M p i = (87rG) -1 / 2 ~ 
1. 

In the ordinary slow roll approximation, the curvature 
power spectrum is given by A^, « / since 
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In the GSR approximation, the curvature power spec- 
trum is instead determined by features in the source func- 
tion through 



lnA^(fc) a G(ln77 min )+ f ^ ^-W (k V )G' (\n V ) 

+ ]n[l + i?(*)], 
where integrating G gives 



(4) 



G(lnr ? ) = -21n/ + -(ln/)'. (5) 
Here the non-linear correction is given by 

Ii(fc) = 4j r™^X(k V )G'(\n V ). (6) 

We take 77, nin = 1 Mpc and rj max = 10 5 Mpc which more 
than covers the range observable to WMAP. The window 
functions 
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define the linear and nonlinear response of the curva- 
ture spectrum to G' respectively. Accuracy of the GSR 
approximation requires the nonlinear response, as quan- 
tified by I 2 , to remain below order unity. We call this 
the GSR condition (see Appendix |B|) . 
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FIG. 1: The first 20 principal components of the GSR source 
G as a function of conformal time to the end of inflation, 
in order of increasing variance from bottom to top. 20 PC 
components suffice to represent inflationary features observ- 
able to WMAP that vary no more rapidly than ~ 1/4 of an 
efold. Here and below, dashed red lines represent power law 
conditions with zero amplitude in the PC components. 
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B. Principal Component Analysis 

The principal components of the WMAP7 Fisher ma- 
trix provide an efficient basis with which to decompose 
the source function ([!]) 



N 



G'(ln?7) = 1 - n s + ^m a 5 a (ln7?) 



(8) 



where the eigenfunctions S a are constructed following 
Ref. [29| by sampling at a rate of 10 per decade in 
77 or equivalently 4.3 per efold of inflation across 77 = 
[1 — 10 ]Mpc. In terms of the width of features in the po- 

1 /2 

tcntial, this limit corresponds to A0 > /4.3. This rate 
is sufficient to capture models that describe the glitches 
in the WMAP7 power spectrum (see ^IVI for a discussion 
of the limitations imposed by the sampling). 

The amplitudes m a then can be incorporated into a 
Markov Chain Monte Carlo (MCMC) likelihood analy- 
sis of the WMAP data. As described in Appendix [5] we 
slightly modify the original approach [29j to improve the 
convergence properties of the MCMC analysis. Since a 
constant G' described by n s is equivalent to tilt in the 
curvature spectrum and G(ln?7 m i n ) is equivalent to a nor- 
malization parameter we replace them with effective pa- 
rameters G' and A c . Specifically G' is an average of G' 
for 30 < 77/Mpc < 400 and A c is the normalization of 
the temperature power spectrum Cj T at the first peak 
I = 220 relative to a fiducial choice that fits the WMAP 7 
data. From these two phenomenological parameters we 
can derive constraints for the tilt n s and curvature power 
spectrum normalization A s (sec Appendix lA)). 

Since a signal-to-noise analysis shows that 20 out of 
the 50 principal components are required for a complete 
representation of the WMAP data at our bandlimit [2^] , 
we choose N = 20 for our analysis. These first 20 prin- 
cipal components are shown in Fig. [1] Note that the 
first 10 components resemble local Fourier modes around 
tj k 10 2 Mpc where the well-constrained first acoustic 
peak gets its power. It is not until components 11-20 
that horizon scale features at low multipole or 10 3 — 10 4 
Mpc are represented. 

We use the MCMC method to determine joint con- 
straints on the 20 PC amplitudes and cosmological pa- 
rameters 

Pfl = {mi, . . . , m 20 , A c , G, r, n b h 2 , Cl c h 2 , 0} . (9) 

Here r is the reionization optical depth, Q(,h 2 is the physi- 
cal baryon density, Sl c h 2 is the physical dark matter den- 
sity and 9 is the angular size of the sound horizon at 
recombination. 

The MCMC algorithm samples the parameter space 
evaluating the likelihood £(x|p) of the data x given each 
proposed parameter set p (e.g. see [3(1 HJ). The poste- 
rior distribution is obtained using Bayes' Theorem, 




FIG. 2: Sensitivity of the nonlinearity parameter 7i imax (see 
Eq. Hip to the amplitude of the first 20 PCs considered indi- 
vidually. This parameter must be less than order unity for the 
GSR approximation to be accurate, and we typically place a 
prior of ii, max < l/\/2- 



where V(p) is the prior probability density. We place 
non-informative tophat priors on all parameters in 
Eq. © • To ensure the validity of the GSR approximation 
we also place a tophat prior on 



(11) 



7>(p|x) = 



£(x|p)7?( P ) 

J de c( x \p)v(p) 



(10) 



h,max = max|ii(fc)|. 



As shown in Appendix [Bl a value of ii jmax = l/\/2 is 
sufficient to ensure accuracy of the GSR approximation. 

Fig. [2] shows the maximal contribution to I\ per unit 
amplitude deviation in each of the first 20 principal com- 
ponents. The higher PCs actually produce a slightly 
smaller response largely because the frequency of the os- 
cillations in Fig. Q] begins to exceed that of the nonlinear 
response function X(krj). Thus a prior of ii. max = l/v2 
actually allows high PC components to reach order unity 
and \G'\ to reach ~ 4 or greater. 

The MCMC algorithm generates random draws from 
the posterior distribution that are fair samples of the like- 
lihood surface. We test convergence of the samples to a 
stationary distribution that approximates the joint pos- 
terior density P(pjx) by applying a conservative Gelman- 
Rubin criterion [32| of R — 1 < 0.01 across four chains. 
We use the code CosmoMC [33[ for the MCMC analysis 

For the WMAP 7 data [3J], we use the optimized ap- 
proximate likelihood from Ref. [29|. In addition, we uti- 
lize data from BICEP and QUAD which include polar- 
ization constraints [HI, [36| . We calculate the CMB power 
spectra with gravitational lensing turned off and the de- 
fault sparse sampling in I (accuracy _boost=l). We cor- 
rect for these approximations in postprocessing by im- 
portance sampling as described in Appendix [5] before 
presenting the results in the next section. The main ef- 
fect is a ~ 0.5cr upwards shift in the fl^h 2 posterior to 
compensate the smoothing effect of lensing. 
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Parameters 


All Data 


CMB Only 




2.241 ± 0.048 


2.233 


2.231 ±0.051 


2.229 


Q. c h 2 


0.1101 ± 0.0040 


0.1098 


0.1110 ±0.0051 


0.1116 


9 


1.0398 ± 0.0022 


1.0397 


1.0394 ± 0.0022 


1.0394 


T 


0.089 ±0.014 


0.086 


0.087 ±0.015 


0.085 


n s , 1 - G' 


0.9669 ± 0.9882 


0.9649 


0.9620 ± 0.0078 


0.9622 


ln[10 1( M s ] 


3.0808 ± 0.0332 


3.0733 


3.0770 ± 0.0338 


3.0746 


Ho 


71.23 ± 1.74 


71.23 


70.70 ± 2.32 


70.40 




0.738 ± 0.0198 


0.739 


0.732 ± 0.027 


0.730 


-21n£ 


8140.06 


7608.39 



TABLE I: Power law (PL) parameter results: means, standard deviations (left subdivision of columns) and maximum likelihood values 
(right subdivision of columns) with CMB data (WMAP7 + BICEP + QUAD) and all data (+UNION2 +H + BBN) in a flat universe. 
Ho and constraints are derived from the other parameters. 



Parameters 


All Data /i max 


= 1/V2 


All Data ii, max 


= 1/2 


CMB Only Ji, max 


= 1/V2 


lOOlW 


2.279 ±0.107 


2.227 


2.282 ±0.107 


2.410 


2.160 ±0.159 


2.110 


tl c h 2 


0.1127 ± 0.0055 


0.1101 


0.1126 ±0.0056 


0.1100 


0.1297 ±0.0142 


0.1338 


8 


1.0411 ± 0.0030 


1.0402 


1.0411 ±0.0030 


1.0417 


1.0395 ± 0.0032 


1.0381 


r 


0.086 ±0.016 


0.096 


0.088 ±0.016 


0.091 


0.082 ±0.016 


0.072 


G' 


0.0122 ± 0.0268 


0.0055 


0.0191 ± 0.0248 


0.0213 


0.0186 ±0.0283 


0.0221 


m[10 10 yl c ] 


0.0032 ±0.0117 


0.0036 


0.0032 ±0.0122 


0.0098 


0.0051 ±0.0122 


0.0056 


7711 


0.0048 ± 0.0073 


0.0060 


0.0025 ± 0.0071 


0.0068 


0.0021 ± 0.0078 


0.0009 


777,2 


0.0152 ±0.0122 


0.0163 


0.0120 ±0.0122 


0.0109 


0.0086 ± 0.0137 


0.0104 


7713 


-0.0120 ±0.0181 


-0.0042 


-0.0140 ±0.0179 


-0.0085 


-0.0151 ±0.0191 


-0.0161 


7714 


0.0427 ± 0.0190 


0.0460 


0.0327 ±0.0171 


0.0481 


0.0455 ±0.0195 


0.0583 


7715 


0.0198 ± 0.0256 


0.0050 


0.0168 ± 0.0249 


0.0486 


0.0165 ±0.0272 


0.0165 


7716 


-0.0156 ± 0.0325 


-0.0089 


-0.0142 ± 0.0328 


-0.0166 


0.0062 ± 0.0377 


-0.0120 


7717 


-0.0061 ± 0.0354 


-0.0015 


-0.0060 ± 0.0333 


-0.0060 


-0.0174 ± 0.0383 


-0.0324 


7718 


0.0278 ± 0.0486 


0.0285 


0.0403 ± 0.0464 


0.0431 


0.0174 ± 0.0505 


-0.0061 


7719 


-0.1239 ± 0.0731 


-0.1436 


-0.0970 ± 0.0670 


-0.1458 


-0.1319 ±0.0770 


-0.1184 


77110 


0.0336 ± 0.0609 


0.0219 


0.0282 ± 0.0602 


0.0462 


0.0150 ± 0.0647 


0.0441 


77111 


0.0759 ± 0.0908 


0.0225 


0.0599 ± 0.0847 


0.0364 


0.0591 ± 0.0966 


0.1339 


77112 


-0.0917 ±0.1027 


-0.1604 


-0.0702 ± 0.0946 


-0.1477 


-0.1100 ±0.1076 


-0.2137 


77113 


-0.0947 ± 0.1129 


-0.1895 


-0.0764 ± 0.1036 


-0.1577 


-0.0506 ±0.1194 


-0.2300 


77114 


0.1116 ±0.1616 


0.2069 


0.0561 ±0.1450 


0.2126 


0.1507 ±0.1714 


0.2103 


77115 


-0.0199 ± 0.2042 


0.0617 


0.0191 ±0.1864 


-0.0091 


-0.0255 ±0.2152 


0.0686 


771l 6 


0.1006 ±0.0975 


0.1318 


0.0837 ± 0.0964 


0.1102 


0.1481 ±0.1043 


0.0772 


77117 


-0.1253 ± 0.2688 


-0.1953 


-0.1094 ±0.2326 


-0.1302 


-0.0575 ± 0.2833 


-0.1376 


77118 


-0.5089 ± 0.2938 


-0.6131 


-0.3322 ± 0.2475 


-0.3798 


-0.4894 ± 0.3083 


-0.6610 


77119 


0.2239 ± 0.3773 


0.2737 


0.1524 ±0.3028 


0.1785 


0.2406 ± 0.3878 


0.5228 


77120 


-0.0742 ± 0.4070 


0.0011 


-0.2472 ± 0.3173 


-0.1789 


-0.1265 ±0.4065 


-0.0113 


71a 


1.0299 ± 0.0671 


1.1296 


1.0075 ±0.0515 


1.0535 


1.0191 ±0.0672 


1.0823 


ln[10 1( M s ] 


3.0387 ± 0.0582 


3.0358 


3.0446 ± 0.0573 


3.0654 


3.0684 ± 0.0626 


3.0726 




71.03 ±2.28 


71.22 


71.08 ±2.28 


73.28 


63.86 ± 5.88 


61.35 


n A 


0.730 ± 0.026 


0.739 


0.731 ± 0.026 


0.750 


0.614 ±0.105 


0.588 


2AlnL 


16.85 


14.26 


17.2 



TABLE II: 20 principal component (PC) parameter results: means, standard deviations (left subdivision of columns) and maximum 
likelihood (ML) values (right subdivision of columns). Fiducial results are for all data and nonlinearity prior 7i ima x = l/\/2, left columns, 
with variations shown in center and right columns. Parameters n s — Q A are derived from the chain parameters. The difference in likelihood 
2A In L is given for the ML values and taken with respect to the corresponding PL maximum likelihood model in Tab. [U 
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FIG. 3: Constraints on the 20 PC amplitudes from the all- 
data analysis with an Ii >max < l/v2 prior. The only signifi- 
cant deviation from the m a = PL expectation (red dashed 
line) is mi = 0.0427 ± 0.0190. The impact of the prior can be 
visualized by taking the maximum amplitude of an individ- 
ual m a that satisfies the prior (solid lines), which implies that 
only mn — m2o are significantly prior limited. The maximum 
likelihood model is shown as starred points. 



In order to ensure that models are compatible with 
a reasonable cosmology we add non-CMB constraints 
from the UNION2 supernovae dataset H0l J _the SHOES 
Hq = (74.2 ±3.6) km/s/Mpc measurement [37| and a big 
bang nucleosynthesis constraint of VL^h 2 = 0.022 ± 0.002. 
These data mainly constrain the energy density compo- 
nents of the universe rather than the inflationary initial 
conditions. We call the combination of CMB and exter- 
nal data the "all data" analysis. We address the impact 
of the -Z^max prior and the non-CMB data in £1111 Bl bclow. 



III. MCMC RESULTS 

In this section, we present the results of the Markov 
Chain Monte Carlo (MCMC) analysis in the principal 
component (PC) space of the GSR source function. We 
discuss the results of our fiducial all data analysis in 
mil Al and address the impact of priors and non-CMB 
data in TOTBl 



A. All Data 

For our fiducial results we use the all-data combina- 
tion of CMB and external data described in the pre- 
vious section. To establish a baseline for the PC re- 
sults we start with the m a = power law (PL) case, 
A^ = A s (k/k p ) ns ~ 1 . Table U gives the mean, standard 
deviation of the posterior probabilities, and the maxi- 
mum likelihood (ML) parameter values for the power law 
model. 
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FIG. 4: Parameter probability distributions from the all-data anal- 
ysis in a flat universe with 7i !ma , x = l/\/2. Dashed lines represent 
the posteriors with approximations for the low i polarization likeli- 
hood and Cg accuracy used to run the MCMC (see Appendix); solid 
lines represent posteriors corrected by importance sampling. Red 
dashed lines represent corrected posteriors for power law models. 



For the PC analysis, we take 20 components and a 
nonlincarity tophat prior of Ji, mnx < \j\Pi (sec < JII Bl) . 
Table [TT] gives the parameter constraints as well as the 
maximum likelihood PC model (left columns). 

The improvement in the ML PC model over the ML 
power law model is 2Aln£ = 17 for 20 extra parame- 
ters and so is not statistically significant in and of itself. 
Of course, specific inflationary models may realize this 
improvement with a smaller set of physical rather than 
phenomenological parameters (see ^IVB[) . and so it is 
interesting to examine more closely the origin of this im- 
provement. 

The main improvement comes from £ < 60 in the TT 
part of the WMA P likelihood with a 2Aln£ = 11.9. We 
shall sec in mV Al that these improvements are largely as- 
sociated with known features in the WMAP temperature 
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power spectrum. 

In terms of the principal components, the improve- 
ments are localized in only a few of the 20 parameters. 
Fig. El plots these m a constraints and ML values. Most 
of the components are consistent with zero at the ~ la 
level. Components myj — 77720 are constrained in part by 
the /i jma x prior not just the data. 

As in the 5 PC analysis of (29| . the single most dis- 
crepant parameter between the PL and PC cases is 7774 
corresponding to a feature centered around r\ ~ 300Mpc 
and resembling a local running of the tilt. Fig. @] shows 
the posterior probability distributions of the parameters. 
An 7774 value as extreme as the power law value of to 4 = 
is disfavored at 98.2% CL compared with 94.8% for 5 
PCs and WMAP7 alone. The increase in significance by 
a fraction of a a arises because of the correlation between 
m4 and the higher principal components. Perhaps more 
importantly freedom in the higher PCs allows large 777,4 
without the need to make large adjustments to the cos- 
mological parameters that would violate non-CMB con- 
straints. On the other hand, one event out of 20 showing 
a 98% exclusion is not that unlikely. 

The poorly constrained a > 10 modes allow large 
amplitude deviations and in fact even marginally pre- 
fer them. This explains why including the higher com- 
ponents can change results on the lower components. 
Large amplitude deviations in the high order compo- 
nents make the modes no longer statistically indepen- 
dent as they would be for infinitesimal deviations. Still 
the correlation remains relatively small. For example 
i?4 Q = Cov(m4, m a )/a mi a ma reaches 0.4 only for one 
mode, 7775, with more typical correlations in the ±0.1—0.2 
range. 

The next most significant deviations are in 7779 (with 
a value of 7779 = disfavored at the 89.6% CL) and mis 
(with a value of mis = disfavored at the 91.8% CL). 
These results are also consistent with the PL null hypoth- 
esis of 777 a = 0, given that there are only 3 events out of 
20 where tests of that model exceed the ~ 90% CL. 

We can get further insight on the origin of these con- 
straints by examining the maximum likelihood (ML) 
models. Fig.[5]show temperature and polarization power 
spectra of the ML PL (red dashed lines) and PC (thick 
solid curve) models respectively. The poorly constrained 
a > 10 modes create fluctuations in the low order multi- 
poles which marginally fit features in the data better such 
as the low quadrupole and glitch at I ~ 20 — 40. These 
large amplitude modes require small amplitude low or- 
der PC variations in order to compensate the broad band 
residual effects they have. This can be seen by decom- 
posing the difference between the ML PL and PC mod- 
els into contributions from the various parameters (see 
Fig. [HI) - Removing the large mio — 77720 components from 
the model not only removes the low £ oscillations but also 
creates broadband deviations, especially at £ < 40, that 
are compensated by a combination of small amplitude 
changes in mi — 7775 and effective tilt G . 



B. Robustness Tests 

In order to test the robustness of the fiducial results 
of the last section, we run separate MCMC chains with 
different choices for the nonlinearity prior and data sets. 

We first examine the impact of our Ii >maK prior by 
reanalyzing the all-data case with Ii >maK = 1/2 instead 
of 1/ \/2 (see Tab. |TT| . The main impact of tightening the 
prior is on m 18 — 777 2 o as is expected from Fig. [3] These 
components mainly affect the low £ multipoles. In spite of 
this fact the prior on ii. max has very little impact on the 
behavior of favored models at low I. In Fig. [Jj we show 
the maximum likelihood model with the stronger ii. max 
prior. Even at low £ the differences are much smaller than 
cosmic variance. In particular the posterior distribution 
of power in the quadrupole moment for models in the 
chain shown in Fig. [8] differ negligibly 

Some of this robustness in the low multipole moments 
is due to the impact of the non-CMB data. Without 
the external data, the quadrupole distribution extends 
to smaller quadrupole moments due to the ability to re- 
duce the integrated Sachs- Wolfe effect by lowering the 
cosmological constant in the absence of constraints on 
the acceleration of the expansion (see Fig. [8]). In this 
case the data may prefer more extreme inflationary mod- 
els that further lower the quadrupole that are excluded 
by our nonlinearity prior on Ii >laax (HI- 

The main impact on parameters of removing the non- 
CMB data is to allow a wider range in f2 c /i 2 (see Table 
HI]) . In contrast to the 5PC analysis [2{|, this wider range 
though has little impact on the PC parameters. In partic- 
ular the higher order PC components allow compensation 
of the effects of 7774 across the acoustic peaks without the 
need to vary fl c h 2 substantially. For similar reasons, we 
expect our flatness prior to have little impact on the PC 
results aside from weakening the constraints on J7a and 
Sl c /7, 2 and small shifts of the location of features in G 
with the angular diameter distance degeneracy. 



IV. APPLICATIONS 

Here we discuss applications of the fiducial 20 PC anal- 
ysis of ^III Al In ^IV Al we place constraints on and de- 
vise tests of slow roll and single field inflation in a model 
independent manner. Alternately, as a complete obser- 
vational basis for efold bandlimitcd models, the PC anal- 
ysis places constraints on any such model that satisfies 
the GSR condition. We use running of the tilt and a step 
in the inflaton potential as example test cases in flIV Bl 



A. Testing Slow Roll and Single Field Inflation 

Bounds on the PC components can be thought of as 
functional constraints on G itself across the observed 
range from WMAP. These in turn limit features in the 
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FIG. 5: The temperature (left) and _B-mode polarization (right) power spectra posterior using the all-data PC constraints and 
a prior of Ji, max = l/s/2. The shaded area encloses the 68% CL region and the upper and lower curves show the upper and 
lower 95% CL limits. The maximum likelihood (ML) model is shown as the thick black central curve, and the power law ML 
model is shown in red dashed lines. The blue points with error bars show the 7-year WMAP measurements. 




FIG. 6: Parameter decomposition of the temperature power 
spectrum difference between the power law (PL) and PC max- 
imum likelihood (ML) models shown in Fig. [5] (bottom panel). 
The curves include cumulative changes in parameters between 
the models starting with the cosmological parameters, adding 
the normalization A c and effective tilt G' , mi . . . ms, etc. un- 
til the full PC ML parameters are utilized. 



FIG. 7: Comparison of the maximum likelihood models of 
the three MCMCs of Tab. [II] the all-data analysis with 
7i,max = l/v2 (black curve), all-data with 7i >max = 1/2 (blue 
curve), and CMB data with Ji imax = l/y/2 (red curve). The 
smallness of the differences indicates robustness of our results 
to the priors and external data sets. 



If the infiaton carries non-canonical kinetic terms then 
inflaton potential V((j)) through the approximate relation the relationship is modified to include variations in the 
[H[ sound speed (27j . 

Since the PC decomposition only represents features in 

V 4> Vit>4> across the observable domain, one should consider the 

G'Qnr,) « 3(^) 2 -2^. (12) 

constraints on the m a s as defining a PC filtered version 
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FIG. 8: The temperature quadrupole power C 2 distribution 
for the all-data analysis with Ji, ma x = l/v2 (black curve), 
all-data with 7i, ma x = 1/2 (blue curve), and CMB data with 
ii.max = 1/V2 (red curve). Without external data to con- 
strain the cosmological constant, the quadrupole can be low- 
ered by reducing the integrated Sachs- Wolfe effect. 

of G": 

20 

G 20 (lnr ? ) = ^m Q 5 a (ln7 ? ). (13) 

0=1 

Any significant deviation from zero of this function would 
indicate a violation of ordinary slow roll. We can extract 
the posterior probability of G 20 by considering its values 
on a continuous set of samples of 77 as derived parameters. 

In Fig. Owe plot both the ML model and the 68% and 
95% posterior bands. Note that G' 2Q = lies within the 
95% CL regime for all 77. These functional constraints 
differ from a full reconstruction of G' in that the PCs 
filter out deviations at 77 < 20 Mpc and rj > 10 4 Mpc as 
well as deviations that are too high frequency to satisfy 
our bandlimit. 

In the well-constrained regime of 30 < 77/Mpc < 400 
constraints are both tight and consistent with G' 2Q = 0. 
Only nearly zero mean high frequency deviations are al- 
lowed in this regime. Nonetheless, the poorly constrained 
toio — 772.20 components allow, but do not strongly prefer, 
large oscillatory features between 10 3 < 77/Mpc < 10 4 . In 
fact G' 2Q = lies noticeably outside the 68% CL bands 
only for the dip and bump between 1000 — 2000 Mpc and 
a bump at 70- 100 Mpc. 

We can associate the most significant features with 
the corresponding effects on the observable power spec- 
tra themselves. Figure [5] shows the 68% and 95% range 
in the power spectra posterior. The 1000 — 2000 Mpc 
feature in fact corresponds to the I = 20 — 40 dip and 
bump in the temperature power spectrum. The 70—100 
Mpc feature corresponds to a glitch at t ~ 600 — 700 

[HOI. While the ^ ~ 104 Mpc 

regime is limited by our 
priors on the amplitude of deviations through ii. ma x we 
have shown that the data do not favor a feature corre- 



sponding to a low quadrupole I — 2 unless acceleration 
constraints are omitted (see £|III Bp . 

Finally, we can examine the posterior distributions of 
the E-mode polarization. These predictions are not sig- 
nificantly constrained by the polarization data sets em- 
ployed. Instead these distributions are limited mainly by 
the common origin of the temperature and polarization 
spectra from single field inflation. These serve as pre- 
dictions for future measurements. For example, the low 
significance features in the temperature power spectrum 
predict corresponding ones in the -E-mode polarization 
which have yet to be measured and can be used to test 
the hypothesis of their inflationary origin at substantially 
higher joint significance .38J. In particular, one expects 
a ~ 26%^ 43 ^ enhancement in the EE power spectrum 

at £ = 39 and a 37%+^T deficit around I = 25. The 

skew distribution in the latter case reflects the difficulty 
in constructing models with low power out of the prin- 
cipal components rather than the data disfavoring such 
models. Models that actually explain the low TT power 
at I = 25 predict low EE power as well. 

Even in the acoustic regime where the polarization 
predictions are tight and do not suggest the presence 
of features, these predictions are of interest. If future 
observations violate them, then not only will slow-roll 
inflation be falsified but all single field inflationary mod- 
els, including those with sound speed variations, as long 
as they satisfy our weak prior constraint on acceptable 
models: the efold bandlimit and small GSR non-linearity 
Ii.max < l/\/2- Such a violation might indicate other de- 
grees of freedom breaking the relationship between the 
temperature and polarization fields, e.g. isocurvature 
modes in multifield inflation or trace amounts of cosmo- 
logical defects. For I < 30 violation could alternately 
indicate a more complicated reionization scenario (38[. 

Currently these bounds and tests apply to the I < 800 
regime measured by WMAP but will soon be extended 
by high resolution ground based experiments and Planck. 

B. Constraining Inflationary Models 

We can also apply the model independent principal 
component analysis to any specific set of models that 
satisfy the GSR condition Zi.max < 1/V2 and bandlimit 
of features no sharper than about 1/4 efold. To place 
constraints on the parameters of a model, one projects 
the source function G' of the model onto the principal 
components 

m a = - f^^Gn^G'Gnr?) (14) 

In 77 max - In 7i rain J Vm . n 77 

as a function of parameters and compares the result to 
the joint posterior probability distributions of the com- 
ponents. Likewise one can construct G' 20 from the result 
and compare it with Fig. [5] 

In fact, the means and covariancc matrix C of the com- 
ponents m a form a simple but useful representation of the 
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FIG. 9: The 20 PC filtered G' posterior from the fiducial all- 
data analysis and Ii, m ax = l/v2 as a prior. The shaded area 
encloses the 68% CL region and the upper and lower curves 
show the upper and lower 95% CL limits. The maximum 
likelihood is shown as the thick black central curve, and the 
power law ML model is shown in red dashed lines. 
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FIG. 10: A model with a linear deviation in G' with slope 
a = —0.026 (and arbitrary offset) is shown as the blue curve. 
The 20 PC filtered source G' 2 o (in black lines) is compared 
with the input linear G' model. 20 PCs captures all of the 
observable information in a. These models are compatible 
with the 68% CL region (shaded) for G'20 from the fiducial all 
data analysis. 



joint PC posteriors. From these, one can construct a x 2 
statistic 



20 

X 2 = [( m a - m a )C~l (m b - rh b )] , (15) 

a,b=l 

or the likelihood L oc exp(— x 2 /2) under a multivariate 
Gaussian approximation to the posteriors. For exam- 
ple the ML PC model gives an improvement of A\ 2 = 
— 15.36 over PL to be compared with — 2A lnL = —16.85. 

As a simple illustration of a concrete model, consider 
a linear deviation in G' 



G'Qnrj) = 1 - n + a In (77/770) 



(16) 



The curvature power spectrum for this model has a local 
tilt of 



d In A 2 , / krio 

= n - 1 + aln ' 



an I\ 

where C = e 7/z -' tE /2 f« 2.895 and 



— (1 — no — a\nkr/o) + 1.67a 



, (17) 



(18) 



For I no — 1| <C 1 and |a| -C 1, the 1\ term contributes 
negligibly and the model gives a linear running of the tilt 

The 20 PC components are a linear function of a given 
explicitly by 



m a (a) 



In 77„ 



In 77 



— 5 a (ln77) ^(77/770) . 

min J rjmi n V 

(19) 

In Fig. [TU] we show an example with a = —0.026 and 
compare the original linear G 1 to the PC filtered G' 2Q - 



The filter introduces features at low and high 77 that are 
not present in the actual source. Note that a Fisher anal- 
ysis of sensitivity to a reveals that most of the signal to 
noise should lie in the 777,4 component [29| which carries 
the most significant deviations from zero in the data. 

The x 2 analysis with all data implies a = —0.039 ± 
0.019. We can compare this result to a direct MCMC 
analysis with a as a parameter constructed from 20 PCs: 
a = — 0.027 ± 0.021. Thus the simple \ 2 approximation 
captures the information on a in the 20 PC posterior to 
- 0.5cr. 

We can further test the completeness of the 20 PC 
decomposition of a by going to 50 PCs. In this case 
a = -0.026 ± 0.023 showing that 20 PCs completely 
describe the observable properties of a. In fact, 5 PCs are 
enough to describe the observable properties of a in this 
case; a direct MCMC analysis gives a = -0.026 ± 0.020. 
Fig. [TT]shows that the full posterior distributions of a for 
these cases are indistinguishable within the errors. We 
also show the simple \ 2 approximation which is shifted 
by ~ 0.5cr as expected. 

The posterior distributions are skewed to negative val- 
ues of a. For example the ML model of the 50 PC 
chain has a = —0.021 to be compared with a mean 
of —0.026. For large negative a, the linear G' model 
no longer matches a running of the tilt due to the I\ 
terms in Eq. (| 18|) . In Fig. [12l we show an example with 
?7 = 0.96 and 770 = 145 Mpc for a = dn s /d\nk = -0.09 
and —0.02. While the a model closely matches constant 
dn s /d\nk for the smaller value, it produces substantially 
less deviations at high and low k. This bias explains 
the difference between constraints on the linear a model 
and running of the tilt found in [2j|. For example, with 
the same data sets and priors running of the tilt gives 
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FIG. 11: Posterior probability distribution of a from a direct 
MCMC analysis constructed from 50 PCs (black/solid curve), 
20 PCs (red curve), and 5 PCs (blue curve). The distribution 
from the x 2 approximation is shown in black/dashed curve. 
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FIG. 12: Initial curvature power spectrum of a model with 
running of the tilt (dn a /dlnk — —0.02, —0.09, solid curves) 
compared to a model with a linear deviation in G' (a = —0.02, 
—0.09, dashed curves). For the —0.02 case, the two models 
are similar whereas for —0.09 the running of the tilt model 
has larger deviations from scale free conditions at low and 
high k. 



dn s /d\nk = -0.018 ± 0.019. Note the ML a = -0.021 
from the 50PC chain is consistent with this constraint. 

Another example is the step potential which has been 
employed to explain the glitches in the power spectrum 
at t w 20 - 40 



V{4>) 



1 



2 j.2 

m , <p 



1 + c tanh 



(20) 



For simplicity, we fix <fi s = 14.668 so that the feature 
appears at the correct position to explain the glitches. 
Although we set the smooth part of the potential to cor- 
respond to an m 2 (p 2 model with m ~ 7.126 x 10~ 6 for 
the projection onto PCs, in the analysis we retain the 
freedom to adjust the amplitude and tilt as usual. This 
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FIG. 13: Constraints on the step potential model parame- 
ters c (height of step) and d width of step. Top panel: the 
X 2 approximation (black curves) compared to the full 20 PC 
posterior (blue curves). Bottom panel: constraints from the 
20 PCs posterior (blue curves) compared to a direct GSR cal- 
culation of the model (black points) . 



leaves us with 2 additional parameters c and d to control 
the amplitude and width of the step. 

The constraints on (c, d) from the y 2 approximation 
are shown in Fig. Q1J] (top panel). Note that the crude 
X 2 analysis correctly picks out the favored parameters 
which can explain the glitches [3fij . The minimum x 2 
model is c = 0.0015, d = 0.026 and is favored over the 
PL m a = (or c = 0) model by Ax 2 = -10.2. Although 
the x 2 analysis assumes that the joint posterior in m a is a 
multivariate Gaussian, it does not make that assumption 
for parameter probabilities. With the distorted shape of 
the confidence region, the 68% contour corresponds to 
Ax 2 = 2.5, 95% contour to 8.6 and 99.7% contour to 
13.3, compared with the more stringent 2.3, 6.2 and 11.6 
obtained for Gaussian distributions in (c, d). Here and 
below we take a prior of d > 0.005 due to our bandlimit 
of 1/4 efold (see below). 

We again compare this with a full analysis of the joint 
20 PC posteriors. As in the case of a, the projection 
onto the two dimensional m a (c,d) space leaves us with 
too few samples in the original 20 PC chain to reliably ex- 
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FIG. 14: The ML model of the step potential from the \ 2 
approximation is shown in blue dashed lines, and the ML 
model from the projection onto 20 PCs [m a (c, d)] is shown 
in black lines. The step potential model captures the main 
feature seen in the fiducial all-data analysis (shaded 68% CL 
area) . 



tract the posterior via importance sampling. We instead 
run a direct MCMC analysis on the 20 PC description 
with m a (c, d). These results arc shown in Fig. [T3]in blue 
points. The maximum likelihood model has c = 0.0016, 
d = 0.025 and is favored over PL m a = (or c = 0) 
by 2AlnL = 9.1. These values are fully consistent with 
the simple % 2 analysis. This improvement is a substan- 
tial fraction of the total of 17 available to the 20 PCs 
from Tab. Ull and is achieved with 3 parameters: c,d and 
implicitly <p s , the location of the step. 

The filtered G' 2Q source for both the ML and minimum 
X model are shown in Fig. [14] and are consistent with 
the posteriors of the fiducial all data analysis. Further- 
more, the x 2 analysis correctly picks out the best fit re- 
gion and qualitatively recovers its distorted shape. The 
main difference is that the confidence region is slightly 
underestimated. 

Finally, we test the completeness of the 20 PC descrip- 
tion of the step model by conducting a separate MCMC 
with the full function G' directly (see Appendix [B] 
Eqs. IB10IIB131 for details). The maximum likelihood 
model has c = 0.0021, d = 0.029 and is favored over 
PL by 2AlnL = 9.5 As shown in Fig. [TJ (bottom 
panel), the main difference is that the models are more 
tightly constrained at d < 0.01. The features in G' span 
less than ~ 1/4 of an efold for these models and con- 
sequently the 20 PC decomposition is not complete. In 
Fig. [15] (top panel) we show a model with d = 9.2 x 10~ 3 
and c = 4.6 x 10~ 4 represented by the full function G' 
(in black lines) compared to its 20 PCs description (in 
blue/dashed curves). The fractional difference between 
these two constructions is shown in the bottom panel. In 
such models, the oscillations in the temperature power 
spectrum continue to higher £, in this case t ~ 100, and 
are not allowed by the data. 
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FIG. 15: Top panel: step potential model with width d = 
9.2 x 1CT 3 and height c = 4.6 x 1CT 4 represented by the full 
source function G' (in black lines) compared to its 20 PC 
description (in blue/dashed lines). Bottom panel: fractional 
difference between the full GSR description and its 20 PC 
decomposition. The oscillations at i ~ 100 are not captured 
by the 20 PCs. 



This example shows that the main limitation of the 20 
PC analysis is that it is too conservative for models with 
high frequency structure in the source: such models tend 
to be in conflict with the data in ways not represented 
by the principal components. 



V. DISCUSSION 

We have conducted a complete study of constraints 
from the WMAP7 data on inflationary features beyond 
the slow roll limit. Using a principal component (PC) 
basis that accommodates order unity features as fine as 
1/10 of a decade across more than 2 decades of the infla- 
tionary expansion, we find no significant deviations from 
slow roll. Although one component shows a deviation 
at the 98% CL, it cannot be considered statistically sig- 
nificant given the 20 components tested. The maximum 
likelihood PC parameters only improves 2AlnL by 17 
for the 20 parameters added. 

On the other hand, specific inflationary models may 
access this improvement with fewer physical parameters. 
Most of the improvement comes from fitting features in 
the temperature power spectrum at multipoles I < 60 
with the known glitch at 20 < I < 40 comprising a large 
fraction. 

From our analysis, we also extract predictions for the 
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corresponding features in the polarization power spec- 
trum that can be used to test their inflationary origin in- 
dependently of a specific choice for the inflaton potential 
(cf. [13]). In particular, one expects a ~ 26% enhance- 
ment in the EE power spectrum at I = 39 and a ~ 37% 
deficit around i — 25 if the temperature features have an 
inflationary origin. Outside of the range of these low I 
features, the predictions are very precise and any viola- 
tion of them in future observations would falsify single 
field inflation independently of the potential. 

Our constraints can also be used to test any single field 
model that satisfies our conditions. Most of the informa- 
tion from the likelihood analysis is distilled in the means 
and covariance of the principal components themselves 
which we make publicly available [5H . Two models illus- 
trate this encapsulation: a linear source model that ap- 
proximates running of the tilt and a step potential model 
that fits the features at I = 20 — 40. A simple x 2 anal- 
ysis approximates the joint parameter posteriors despite 
its highly non-Gaussian form for the step parameters. 
This procedure greatly simplifies the testing of inflation- 
ary models with features in that parameter constraints 
on any model that satisfies our conditions can be simply 
approximated without a case-by-casc likelihood analysis. 
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Wyman for useful conversations. This work was sup- 
ported by the KICP under NSF contract PHY-0114422. 
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Appendix A: MCMC Optimization 
1. Parameterization 

We seek to define amplitude and tilt parameters for the 
MCMC that are nearly orthogonal to the PC parameters 
in order to improve the convergence properties of the 
MCMC chains. 

A constant G is equivalent to tilt n s and hence PC 
components that have long positive or negative definite 
stretches become degenerate with tilt and cause problems 
for MCMC convergence. Instead of a constant tilt, we 
define a new chain parameter to be the average of G 
across a narrower range that is better associated with 
the obscrvablcs 



G 



1 



In 772 - hi 771 J Vl 7} 



m dn 

— G' , (Al) 



where specifically, we choose r\\ = 30 Mpc and 772 = 400 
Mpc to roughly minimize the variance of G in the chain 
(see Fig. |U). 

Next, we replace the normalization parameter 
G(ln7fa lin ) with 

A s =\nA 2 K (k p ), (A2) 
where in practice we choose k p = 0.05 Mpc -1 . 



The effective tilt and normalization parameters bring 
the model of the power spectrum from Eq. @ to 



lnA^ = In 



N 



In 



In 



+ ^m a \Wa(k)-Wa(kp)] 
a=l 



where 



W a (k) = P"" ^W{kri)(S a {\nri)-S a ), 

«/?7 min V 

X a (k) = ^lx(k V )(S a (ln V )-S a ), (A4) 



and 



S a = 



dn 

In ?y 2 - In 771 J m 77 



S a . 



(A5) 



Note that we can recover the tilt n s , equivalent to the 
average of G across the whole range ?7 m i n to 77 ma x, as 



20 



(l-G)+J2m a S a , 



(A6) 



and keep it as a derived parameter in the chain 

Given the oscillatory nature of the fc-space response 
to the PC eigenfunctions through W a and X a and the 
geometric projection from k to angular multipole nor- 
malization at a given k does not correspond simply to 
normalization at a given t. Since the observations best 
constrain the amplitude of the temperature power spec- 
trum near the first acoustic peak at I ~ 220 it is advan- 
tageous to use an ^-space normalization in the MCMC 
and then transform back to A s . 

Let us define a phenomenological parameter A c which 
renormalizes the angular power spectra as 



a 



XY 



„lnA c °220 riXY 

fiTT ° f • 
°220 



(A7) 



Thus if A r = 0, C?t 



Here C|2o fid is the temperature power spectrum at the 
first peak of a fiducial model that fits the WMAP7 data. 
0.747/jK 2 . 
regardless of the PC parameters. 
We can recover constraints on the fc-space normaliza- 
tion by considering A s as a derived parameter. If we 
compute the original Cf Y with the A s = A^ d of the 
fiducial model, then the true A s is given by 



We use C™ 

°220 



In A s = In A c 



ln(C 2 T 2 r/^)+ln^ d . 



(A8) 



In summary, we replace the parameters n s and G(hi77 m i n ) 
with G and A c in order to reduce parameter degeneracies 
that would otherwise inhibit chain convergence. 
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2. Likelihood corrections 

To speed up the calculation of the WMAP and other 
CMB likelihoods we employ three approximations when 
running the chains. Firstly, we use a fitting function for 
the low I WMAP7 polarization likelihood as described 
in [2^]. Secondly, we calculate the CMB power spectra 
with gravitational lensing artificially turned off. Thirdly, 
we use the default £-space sampling of CAMB that is 
designed for smooth underlying power spectra. Each of 
these approximations produce small errors in the likeli- 
hood evaluation that we can correct by importance sam- 
pling the chain. 

The advantage of correcting these approximations in 
a postprocessing step is twofold. The chains may be 
thinned due to the high correlation between samples 
in the chain. Secondly, postprocessing elements of the 
thinned chains is embarrassingly parallel unlike the run- 
ning of the original chain. 

In practice, when we satisfy our convergence criterion 
described in the main text, we thin the chains by a factor 
of half of the correlation length. We have tested that with 
such thinning we reproduce the posteriors of the original 
chains. Next we compute the CMB power spectra of 
the thinned chains with lensing turned on and a higher 
£-space sampling (CAMB "accuracy boost" 2). We use 
these high accuracy power spectra to correct the chain 
multiplicity for the change in the likelihood. 

In Fig. [16] we show as an example the posteri- 
ors coming directly from power law (PL) chains (in 
blue/solid curves), the chains with all corrections (in 
blue/dashed lines) and finally all corrections but lensing 
(in red/dashed lines). These should be compared with 
results from a separate chain run with all the corrections 
turned on from the start (in black/solid lines). Impor- 
tance sampling accurately models the impact of the small 
corrections for all parameters. The leading correction is 
on flth 2 from lensing. 

In Fig. [T71 we show the impact of the corrections on the 
PC chain using as an example with the largest cor- 
rection. The correction on PC parameters is extremely 
small and again dominated by lensing. 



Appendix B: GSR Accuracy 

To test the accuracy of the GSR approximation in the 
PC space, we need to consider the inverse problem: con- 
struct an inflationary model that matches a desired G 
for which we can solve exactly for the curvature power 
spectrum. 

In the forward direction, given an inflationary model 
we can compute the exact curvature spectrum by first 
evaluating the background behavior of the model through 

rfM, s£ = l<?+(f) 2 , (Bi) 
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FIG. 16: Power law parameter posteriors from the approxima- 
tions used to run the MCMC chain (in blue/solid curve), from 
an independent MCMC with no approximation (in black/solid 
curve), from the approximate chain with importance sampling 
correction (in blue/dashed curve), and from the approximate 
chain without lensing correction (in red/dashed curve). 




FIG. 17: The mig probability distributions from the approxima- 
tions used to run the MCMC with all data and 7i jma x = l/v2 
(in black lines) , from the approximate chain with importance sam- 
pling correction (in blue/dashed lines), and from the approximate 
chain without lensing correction (in red/dashed lines), mis has 
the largest correction of the PC amplitudes which is still -C la and 
dominated by the lensing correction. 



and then solving the equation 

^J y= ^' (B2) 

where x = krj, subject to the usual Bunch-Davics initial 
conditions. The curvature power spectrum is then given 
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by 



which has the general solution 



A 2 



lim x —r- 

x-)-0 f 2 



(B3) 
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FIG. 18: Fractional difference in temperature power spec- 
tra between GSR and the exact inflationary solution for 
the maximum likelihood model from the all-data analysis 
(in black lines) as well as a model that saturates the prior 
-Ti.max = l/v2 from the chain (in blue lines). For reference, 
the ML model has Ji, m ax = 0.66. 
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FIG. 19: Likelihood difference between the GSR solution and 
the full inflationary calculation of a series of step potential 
models as a function of 7i,max- Models were chosen from the 
full GSR chain to be the maximum likelihood in a series of 
bins in step amplitude c. The maximal error is small below 
-fi.max = l/v2 (blue dashed line), the prior in the hducial 
all-sky analysis. 

Therefore to test the GSR approximation we hrst need 
to determine the function g that matches a given G"(ln 77). 
Transforming variables to r = f'/f, we obtain from 
Eq. JEQ) 



3r 



(B4) 



Let us choose the integration constant 

c=4 r^c, 



2 ./„ nin 77 



(B5) 



(B6) 



and assume G' vanishes outside this range. We then get 



r = -2 n 



for 77 > 77 min and 



r = -- v 



3„ 3 r-d^-s^ 



(B7) 



(B8) 



for 77 < rjmin- With this numerical solution we construct 
g as 



(B9) 



This suffices to specify the source for y in Eq. (|B2[) . Fi- 
nally, to get the curvature power spectrum we need / at 
some 77u m <C fc"^. However since this quantity is inde- 
pendent of k, it is absorbed into our normalization defi- 
nition. In Fig. [TBI we take parameters from the all-data 
chain and use this technique to calculate the temperature 
power spectra of matching inflationary models exactly. 
Even for the model that saturates the ii, m ax = l/v2 
prior, the WMAP likelihood difference between the ex- 
act and GSR calculation is |2AlnL| = 0.4. 

Using the step model chain from ijlVBl we can explore 
the accuracy of the GSR approximation as a function of 
^i,max independently of the prior taken in the all data 
analysis. Specifically, we take a model from the chain 
that defines G' and construct the matching full inflation- 
ary model as above. 

Recall that to construct G', we solve for the back- 
ground evolution of <j> in the step potential of Eq. (|20p . 
This specifies the m 2 <p 2 model source through Eq. ((TJ), 
which we call G' m . To allow for a retilting of the spec- 
trum, we add an extra constant parameter G' p to the 
model source to form the full source 



G' (In 77; c,d,G' p ) = G^(ln77;c,d) + G' p . 



(BIO) 



The GSR approximation then tells us that the curvature 
spectrum is given by 



In A 



■R 



In 



-G' 



I m {k)~I m (k p ) (Bll) 



with 
I m {k) 



drj 



W(k v )G' m 



(B12) 



In 



drj 



X(k V )G' m 
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In Fig. 1191 wc compare the impact of taking this power 
spectrum to a full inflationary calculation with match- 
ing source ()B10j) on the WMAP likelihood. For the full 
calculation of Eq. (|B3[) . we take A^(fc p ) = A s to define 
the normalization /. Since Ix, max increases monotoni- 
cally with c, we show models with maximum likelihood 
parameters in uniform bins of c. Note that the maxi- 
mal error increases with ii. max but does not exceed order 
unity at h, max < 1/^/2. 



For reference, to compute a matching 20 PC repre- 
sentation as in Fig. [TS] we take the amplitudes of the 
principal components from Eq. (fl"4")) and use Eq. (|A1[) to 
define 

G 1 = G' + f 2 *L G' m (lnrj) . (B13) 

p In j] 2 - In 77i J m J] 

and keep the normalization A s fixed. 
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